The critical temperature of a trapped, weakly interacting Bose gas 
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We report on measurements of the critical temperature of a harmonically trapped, weakly inter- 
acting Bose gas as a function of atom number. Our results exclude ideal-gas behavior by more than 
two standard deviations, and agree quantitatively with mean-field theory. At our level of sensitivity, 
we find no additional shift due to critical fluctuations. In the course of this measurement, the onset 
of hydrodynamic expansion in the thermal component has been observed. Our thermometry method 
takes this feature into account. 
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Degenerate atomic Bose gases provide an ideal test- 
ing ground for the theory of quantum fluids. First, their 
diluteness makes possible first-principles theoretical ap- 
proaches 1]. Second, thanks to the powerful experi- 
mental techniques of atomic physics, static and dynamic 
properties can be studied quantitatively through a wide 
range of temperature and densities. Furthermore, the in- 
homogeneity induced by the external trapping potential 
leads to entirely new behavior, when compared to bulk 
quantum fluids. 

Atomic interactions have previously been found to af- 
fect deeply the dynamical behavior of trapped Bose gases 
at flnite temperatures By contrast, the influence 

of interactions on thermodynamics is less pronounced , 
and has been less studied experimentally. Pioneering 
work on thermodynamics concentrated essentially on 
the ground state occupation, and the role of interactions 
was somewhat hidden by finite size effects 0. Though 
several such measurements have been reported 0, , to 
our knowledge a decisive test of the role of interactions 
is still lacking. In this Letter, we focus on the critical 
temperature Tc of a harmonically trapped ^^Rb Bose gas 
to demonstrate the infiuence of interactions on the ther- 
modynamics. We study the behavior of Tc as a function 
of the number of atoms at the transition, for a fixed trap- 
ping geometry. We find a deviation from ideal-gas behav- 
ior, towards lower critical temperatures, whose significa- 
tion will be discussed below. In the course of this study, 
we have observed that collisions induce an anisotropy in 
the free expansion of the cloud even far from the hydro- 
dynamic regime 0, 1^1 ■ We correct for this effect in our 
temperature measurement. 

For an ideal Bose gas in a harmonic trap, the critical 
temperature is 



The first term on the right hand side is the transition 



temperature in the thermodynamic limit, and the sec- 
ond represents finite-size corrections. Here N is the total 
atom number, and lOz are the trapping frequencies, 
Lu = uj'^j^'^Lul^^ is their geometrical average, and ^ is the 
Riemann zeta function. 

As stated earlier, the main goal of this paper is to 
probe the role of two-body repulsive interactions on Tc- 
Corrections to the ideal gas formula depend on the ra- 
tio between the s-wave scattering length a and Ao — 
h I y^MkBT^, the de Broglie wavelength at the transition 
(Ao = (27r)i/2 cr[7 V/C(3)] ~^/^ for a harmonically trapped 
gas, where a = ^JhjMuj is the mean ground state width). 
In a trapped gas, the dominant effect of interactions can 
be understood using a simple mean-field picture in- 
teractions lower the density in the center of the trap ri(0), 
and accordingly decrease the temperature T that meets 
Einstein's criterion n(0)Ai^ = C(3/2). The magnitude of 
this reduction has been calculated to leading order ^J, 
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_a « _1.326-iVl/^ (2) 
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where oi ~ 3.426 this work, the finite-size cor- 

rection in Eq. Q changes by at most 2%, whereas 
the interactive shift (O can be as high as 10%. The mea- 
surements presented below are in quantitative agreement 
with the prediction of Eq. 

In addition, as dicussed in [1 [H [H [H , crit- 
ical fluctuations that develop in the system near Tc are 
expected to favor the formation of the condensate and 
thus to increase Tc. In the case of a uniform Bose gas 
lUfcilS. JJjj this is the leading effect, because the critical 
temperature is not affected at the mean-field level. The 
correction i5Tc/T° = -|-cia/Ao, with c\ « 1.3 Q, can be 
traced back to density fluctuations with wavelength much 
larger than the correlation radius r^ ^ Ap/a jl3l |. This 
upwards trend, which has been observed experimentall 



:aliy 



*UMR 8501 du CNRS 



in a dilute sample of He adsorbed in a porous glass 
is quite sensitive to the presence of an external poten- 
tial 1^, ^1 . In the harmonically trapped case of interest 
here, the contribution of long wavelength excitations to 
the shift in Tc scales as a higher power of a/Ao, making 
it negligible when compared to the "compressional" shift 
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given by Eq. ((SJ- The quantitative agreement we find 
with the mean-field result can be considered evidence of 
this effect, and highlights the important role played by 
the trapping potential. 

Our experimental setup to reach Bose-Einstein conden- 
sation in the \F = l;m-p = —1) liyperfine ground state 
of ^^Rb is similar to that used in The trapping fre- 
quencies are u±/2n = 413(5) Hz and uJz/2tt = 8.69(2) Hz 
in the present work. To reduce non-equilibrium shape os- 
cillations that occur in such anisotropic traps upon con- 
densation H,0|, the last part of the evaporation ramp is 
considerably slowed down (to a ramp speed of 200 kHz/s) 
and followed by a 1 s hold time in the presence of a radio- 
frequency shield. We ensure good reproducibility of the 
evaporation ramp in the following way. We monitor reg- 
ularly (typically every four cycles) the radio-frequency vq 
that empties the trap. This allows us to detect slow drifts 
of the trap bottom, and to adjust in real time the final 
evaporation radio-frequency v^f to follow them. In this 
way, the "trap depth" Vj-f — I'o, is kept constant within 
±2 kHz. Since we measure rj = h{i'ii — VQ)/k^T w 11 in 
this final evaporation stage, we estimate the temperature 
stability to be ±10 nK. 

We infer the properties of the clouds by absorption 
imaging. After rapid switch-off of the trap (1/e cut-off 
time of about 50/Lts), a 22.3 ms free expansion, and a 
repumping pulse, we probe the ultra-cold cloud on res- 
onance with the |F = 2) ^ = 3) transition 
The images are analyzed using a standard procedure, de- 
scribed for instance in For an ideal thermal cloud 
above the transition point, the evolution of the density 
in time of flight is related to the initial density profile 
by simple scaling relations, so that the column density 
(integrated along the probe line-of-sight, almost perpen- 
dicular to the long axis of the trap) is 

nM = -th(0)32{exp(^-^-^)|(3) 

where g2{u) = Sj>i'""'7i^5 ^'^'^ ^ ^'^'^ ^ the coor- 
dinates along the tight and shallow trapping axes, re- 
spectively. For mixed clouds containing a normal and a 
(small) condensed component, we assume, as usual for a 
condensate in the Thomas-Fermi regime m, that one can 
describe the bimodal distribution by an inverted parabola 
on top of an ideal, quantum-saturated thermal distribu- 
tion (Eq. with /i = 0). The condensed number iVo is 
then deduced from integration of the Thomas- Fermi fit, 
and the total atom number TV from integration over the 
entire image. We estimate that condensed fractions as 
low as 1% can be reliably detected by the fitting routine. 
Absolute accuracy on the value of N and relies on the 
precise knowledge of the absorption cross-section of the 
probe laser, which depends on its polarization and the 
local magnetic field. This cross-section is calibrated by 
fitting the radial sizes of condensates with no discernible 



thermal fraction to the Thomas- Fermi law Rq rx Nq' , as 
explained in 0,[1E1E3- We find a reduction of 4.00(14) 
compared to the reference value ctq — 3X1/2tt [i7| . 

We will now discuss the more complex issue of ther- 
mometry in some detail. The temperature is usually in- 
ferred from the sizes of the thermal cloud after a time 
of flight t, assuming a purely ballistic expansion with 
isotropic mean velocity, vg = ^Jk-^T jM , as appropri- 
ate for an ideal gas. We show in Fig. ^ that the ob- 
served aspect ratio of non-condensed clouds, in a wide 
range of temperatures and atom numbers (correspond- 
ing to 1 < T/Tc < 1.8), is actually larger than the value 
(0.773) expected for an ideal gas and Uzt « 1.23 that 
corresponds to our parameters (dotted line in Fig. , in 
contradiction with the assumption of an isotropic veloc- 
ity distribution. 

In a very elongated trap, this could be explained by 
two distinct coUisional effects. First, the initial mean- 
fleld energy of the non-degenerate cloud converts almost 
completely into radial kinetic energy during time-of-fight 
[iiij l , as for an elongated condensate 0, 0| ■ The mag- 
nitude of this effect is controlled by the ratio x of the 
mean-field energy to the temperature. In our case, the 
parameter x does not exceed 0.02, too low to explain 
the observed anisotropy (dashed line in Fig.^ calculated 
along the lines of [20|). 

Second, as studied theoretically in [3. 121I and ob- 
served in Bose and Fermi gases j23l| anisotropic ex- 
pansion occurs for a cloud in the hydrodynamic regime, 
i.e. when the mean free path at equilibrium is smaller 
than the dimensions of the sample. In our very elongated 
cloud, the mean free path is typically smaller than the 
axial length, but much larger than the radial size. Hy- 
drodynamic axial motion of the thermal particles results 
in energy transfer from the axial to the radial degrees of 
freedom. For weak deviations from ballistic expansion, 
this coUisional dynamic creates a velocity imbalance pro- 
portional to 7coii, the equilibrium collision rate, in agree- 
ment with the trend observed in Fig. 

In ji^, a set of scaling equations was derived to in- 
vestigate how collisions affect the expansion of a non- 
condensed cloud. Numerical solution of these equations, 
that also include the weak mean field effect, agrees well 
with our data (solid line in Fig.^. The calculation makes 
use of the results of 24] for the collision rate of a non- 
condensed, almost ideal Bose gas, in general larger (by 
as much as 70% close to Tc) than the classical collision 
rate with the same N and T 25] . In view of the satisfac- 
tory agreement of our data with the scaling theory, we 
conclude that the observed anisotropy is a signature of 
the onset of hydrodynamic expansion. 

In the regime 7coii ^ where the anisotropy is 
weak and increases linearly with 7coii, kinetic energy 
conservation suggests that mean square expansion ve- 
locities take the form {v'^)/vq « 1 -|- /97co1i/2w_l, and 
{Vz)/vq « 1 — /?7cou/^^±, where P depends in general 
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on uj±,ujz,t (fixed for the measurement presented here). 
These simple forms are confirmed by the numerical solu- 
tion described above. Provided the expansion velocities 
along both axes are measured, they allow to infer the ini- 
tial mean square velocity vq and temperature T, indepen- 
dently of the coefficient /3 |2^. Were this correction not 
applied, a systematic 10-15% discrepancy between the 
axial and radial temperature would remain. Note that 
mean-field effects are not corrected by this procedure. 
As stated above, they change the expansion energy by 
2% at most. This error is not significant when compared 
to calibration uncertainties, that limit the accuracy of T 
to about 5% 113 ■ 

Having identified an appropriate thermometric tech- 
nique, we turn to the measurement of the critical tem- 
perature as a function of atom number. Data were taken 
in a narrow range around T^. From the two-component 
fit, we extract the number of condensed atoms, the tem- 
perature, and the total atom number as a function of the 
trap depth, as shown in Fig. b and c, respectively. 
The trap depth at which the transition point is reached, 
(i^rf — 1^0)0, is taken to be the point at which a linear 
fit to the condensed number data crosses zero (a linear 
approach towards Tc is consistent with the simulations 
reported in |28l|V The temperature and total number are 
also fitted assuming a linear dependency on i/j-f , and from 
the value {vii — i'q)c, we extract the critical atom number 
Nc and critical temperature Tc. 

In Fig. 121 we have plotted Tc as a function of Nc, mea- 
sured in ten independent data sets. The ideal gas value 
jndeai (j^shed line) lies two standard deviations above 
our data. Including the mean-field correction ^ yields 
a much better agreement (solid line), that we quantify in 
the following way. We assume that the interactive shift 
in Tc can be written as STc/T^ = aN^^^, with a free co- 
efficient a. A fit to the data yields a = -0.009(l)t^:g^i, 
whereas Eq. ^ predicts a « —0.007 for a scattering 
length a w 5.31 nm [i^ and a « 1.00 ^m. The first 
uncertainty quoted is statistical, while the upper and 
lower bounds reflect calibration and analysis uncertain- 
ties. The shaded area in Fig. |31 dehneates the resulting 
la confidence interval compatible with the experimental 
results. 

The data shown in Fig. |21 reasonably exclude any ad- 
ditional shift of the same order of magnitude as the com- 
pressional effect given by Eq. jSJ- In particular, if the 
(po sitive) critical shift in Tc predicted in the uniform case 
[131 were directly scalable to the trapped one, one would 
expect an overall a « —0.004, a value not consistent with 
our findings within the estimated accuracy. This obser- 
vation is in line with recent theoretical studies of crit- 
ical density fluctuations in a harmonically trapped gas 
0, Q|, which point out that, instead of being delocal- 
ized over the entire system as in the homogeneous case, 
critical fluctuations in the trapped gas are confined to a 
small region around the trap center. This reduces cor- 



rections to the equation of state by a factor ^ (a/Ao)"^, 
corresponding to the ratio of the volume of the fluctu- 
ation region to the volume of the thermal cloud. The 
critical temperature, being fixed by the equation of state 
of the whole cloud, thus depends only weakly on criti- 
cal phenomena: corrections to mean field behavior [s^l 
enter only to second order in a/Ao 0, 0L For our ex- 
perimental parameters, we calculate from |2( an upwards 
correction to Tc smaller than 1%, below the sensitivity of 
the measurement. 

In summary, we have measured the critical tempera- 
ture of a trapped, weakly-interacting ^^Rb Bose-Einstein 
gas. Our results exclude ideal gas behavior by two stan- 
dard deviations, and we find satisfactory agreement with 
mean-field theory. We find no evidence for critical behav- 
ior close to Tc within our experimental sensitivity, in line 
with recent theoretical estimates that predict an increase 
of Tc due to critical fluctuations significantly smaller in 
the trapped case than in the uniform case |^. We have 
also observed hydrodynamic behavior in the expansion 
of the thermal cloud, and shown how to correct for it in 
the thermometry procedure. We note to conclude that 
measuring corrections to Tc beyond the mean-fleld for 
our typical experimental parameters would require ther- 
mometry with an accuracy of 1% or better. A more di- 
rect route to investigate such effects might be to measure 
the critical density near the center of the trap, directly 
sensitive to the presence of critical fluctuations. Alterna- 
tively, these many-body effects could be enhanced in the 
vicinity of a Feshbach resonance |^ . 
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FIG. 1: Onset of hydrodynamic expansion for trapped clouds 
above threshold. Measured aspect ratios after expansion 
(filled circles, with statistical error bars) are plotted versus the 
collision rate at equilibrium 7coii- The experimental results 
are compared against several hypotheses: a ballistic expan- 
sion (dotted line) ; a mean-field dominated expansion (dashed 
line); and a coUisional expansion for a non-condensed Bose 
gas (solid line). 
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FIG. 2: Procedure to locate the transition point. We plot 
the condensed number (a), temperature (b), and total atom 
number (c) as a function of the trap depth, fixed by the final 
rf frequency Vit and the trap bottom uq. The transition point, 
shown as a hollow circle (with statistical error bars) , is found 
from a fit [dotted curve in (a)], and reported in (b) and (c) 
to find Tc and iVc. 
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FIG. 3: Critical temperature as a function of atom number 
at the transition. The experimental points (circles) are lower 
than the ideal gas law Eq. (0 (dashed) by two standard devi- 
ations. The shaded area is the range of acceptable fits taking 
statistical and systematic errors into account. Our results 
are consistent with the shift due to the compressional effect 
given by Eq. indicated by the solid line. The filled circle 
represents the data of Fig. |21 



